NUCLEON-NUCLEON SCATTERING IN 
THE l/N c EXPANSION 



SILAS R. BEANE 

Department of Physics, University of Washington, 
Seattle, WA 98195 
E-mail: sbeane@phys.washington.edu 

The nucleoli- nucleoli 3 Si — 3 D\ coupled-channel problem is solved analytically to 
leading order in a joint expansion in the quark masses and in 1/N C . An approximate 
expression is derived for the 3 Si scattering length in the large-7V c limit, and the 
large- N c behavior of the deuteron is discussed. 

1 Introduction 

When the number of colors, N c , of the QCD gauge group is taken large, QCD 
simplifies. One may thenJiope to learn about the real world with N c = 3, 
in an expansion in l/A^EJfl. The relevance of the large- N c expansion for nu- 
clear physics was first discussed in the seminal paper by Witten, where it was 
argued that the nucleon-nucleon (NN) potential grows as N c in the large- N c 
limit, and Hartree theory becomes exact el. In a recent revival of interest, it 
was shown that large- N c nuclear interactions are spin-flavor SU(4) symmetric, 
with Wigner's phenomenologicalhy successful supermultiplet symmetry follow- 
ing as an accidental symmetry Subsequently, the large-iV c scaling of the 
NN potential was analyzed in a general way and large-iy c expectations were 
shown to explain generic features of the NN interaction □. Given the levels of 
complexity that must be unravelled in going from QCD to a theory of nuclear 
forces, it is quite remarkable that large- N c counting rules work so well. Hence 
large- N c methods offer an important link between QCD and nuclear physics. 
The establishment of such links is, of course, an important goal for nuclear 
physics. In this paper, basic aspects of NN scattering in the large- N c approxi- 
mation will be reviewed and the 3 5i — 3 Di NN coupled-channel system will be 
analyzed in the large- N c limit. I will also say some words about the deuteron 
and the bound-state spectrum in the large- N c limit. 

Assuming confinement, in the large- N c limit the QCD chiral symmetry is 
spontaneously brokenH. With two light flavors the chiral symmetry breaking 
pattern, U{2) (g> U(2) — > U(2), gives rise to four Goldstone bosons, the pions 
and r\. While the chiral symmetry breaking scale, A x , and the meson masses 
scale as N°, the baryon masses scale as N c . In the large- N c limit the nucleons 
and As are degenerate and fall into an infinite-dimensional representation of 
a contracted SU(4) spin- flavor symmetry In this paper I will assume that rj 
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and the As can be integrated out of the low-energy effective field theory_.(EFT) 
of large-iV c QCD in a manner consistent with the large- N c expansion!!!. This 
is, of course, a highly non-trivial assumption. Consider the xesulting EFT 
involving nucleons and pions below A x (for reviews, see Ref. 0). This EFT 
consists of nucleons emitting and absorbing pions with an interaction highly 
constrained by chiral symmetry. The most general Lagrange density consistent 
with QCD symmetries is 



C = \id t + V 2 /(2M)] N- - [C${N^Nf + C^(N^aN) 2 ] 

+ Tr c^E^E + uj^ TrM q (S f + E) + g A N*A ■ aN + . . . (1) 

where A is the axial-vector current and the ellipses represent operators with 
more derivatives, more powers of the quark mass matrix, and more nucleon 
fields. In the large- N c limit, F„ oc \fNc, Qa oc N c and oj oc AT?B. The four- 
nucleon operators scale as Cq oc 1/N c and Cq oc N c u. When considering the 
interactions of more than one nucleon, the ordering of the operators in Eq. (Q) 
by way of a power-counting scheme is nontrivial. This is especially so in the 
3 5i — 3 Z?i coupled-channel of NN scattering in an EFT with pions, which I will 
focus on in this paper. The presence of the deuteron, a near-threshold bound 
state, and the highly-singular tensor force imply complicated nonperturbative 
physics. 



2 NN Scattering at Large- N c 

Defining a sensible large- N c limit for NN scattering is nontrivial because the 
nucleon mass grows as N c . There are several limits we may consider. The 
momentum-space Schrodinger equation is 

+ ^ >=E ^> ' ( 2 ) 

where M is the nucleon mass. At fixed velocity, the momentum transfer scales 
as N c . Since the potential also scales as iV c , an overall factor of N c scales out 
of the Hamiltonian and Hartree theory becomes exact Evidently there is 
no EFT description in this limit as characteristic momenta quickly overcome 
A x and a description in terms of quarks and gluons becomes appropriate. At 
fixed momentum transfer, which is the limit relevant to EFTQ, the nucleon 

a In principle, the momentum transfer can scale as an inverse power of N c in the EFT, but 
then one is constrained to threshold. 
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kinetic energy is suppressed by l/N c while the potential energy grows_as N c . 
In this limit, the nucleons are always below threshold for scattering □. The 
static nucleons settle at-.the bottom of the potential well and nuclear matter 
forms a classical crystal!! The absence of a suitable large- N c limit to describe 
scattering at fixed momentum transfer poses something of a dilemma. One 
way out is to focus on the large- N c scaling properties of the NN potential and 
then to compare the resulting predictions with modern phenomenological NN 
potentials that fit phase shifts to high precision trta. This has proved highly 
successful. The question of whether the fixed momentum^ transfer limit exists 
has been called into question in interesting recent workE9'U which points out 
that there would appear to be violations of the counting rules for the potential 
arising from diagrams with multiple meson exchanges. This is an important 
unsolved puzzle. 



3 The (Un)Coupled Channel Problem at Large-iV c 

The most economical power-counting scheme for the 3 Si — 3 -Di coupled channels 
is an expansion about the chiral limit E3. The pion-exchange part of the NN 
force in the chiral limit is purely tensor, 

V T (r) = -If + 0(ml) , (3) 

where 

~ - 9 " A (4) 
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At leading order in the EFT expansion there is a local four-nucleon operator 
containing the effect of non-pionic short-distance physics which renormalizes 
the tensor force at short distances. The coefficient of this operator is Co = 
Co + Co = c o + 0(1/N c )a. I choose to regulate this coordinate-space delta- 
function using a square well, 

C S^(r)^ 3Co ^- r) =V e(R-r) , (5) 

with strength Vq and width R. The S'-wave and -D-wave components of the 
.7=1 system, u(r) and w(r) respectively, are coupled by the tensor potential, 
Vt, defined in Eq. (0). The long-distance potential outside the square well is 

Vr(r)-( ° -2V2 MV T (r)\ 

LV> \-2V2MV T (r) 2MV T (r) -Q/r 2 ) ' K > 
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while the short-distance potential inside the square well is 



Dehning ^ to be 



u)(r) 

the regulated Schrodinger equation takes the compact form 

*"(r) + (fc 2 + V L (r)6(r - R) + V s (r)0(R - r)) * = 0. (9) 

One may ask whether there is a well-defined large- iV c limit. In Eq. (||), MVr 
and MVo scale as iV 2 , while the angular- momentum barrier scales as N®. In 
the EFT we treat the momentum transfer as fixed (-/V c -independent). The total 
center-of-mass-energy is then suppressed by 1 / compared to the interaction 
energy. Here we will focus on this rather peculiar large- N c limit. At leading 
order in the 1/N C expansion, the coupled channel problem can be diagonalized 
via the similarity transformation ^ — * S~ x ^f where 

s=(f-^). do) 

Clearly 3 Si — 3 Di mixing is suppressed in the large- N c limit. The general solu- 
tion of Eq. (|) in the large- N c limit is a linear combination of Bessel functions. 
Keeping the leading contributions I find 



u(r) = V2A r 3 / 4 cos [ 2 + <*o ) - 4= ^ cosh J^M. + S 



V2 



w(r) = A r 3 / 4 cos + ^ + B r 3 / 4 cosh + 6 ) , (11) 

where (f>o, Sq, A and B are determined from boundary conditions. Note that 
the repulsive hyperbolic solutions grow exponentially with N c . Since these 
solutions cannot match to a power law in N c , we find B = in the large- N c 
limit. We then have 



t(r) - v/2 A r 3 / 4 cos ( 2 J + O J = y/2 w(r), (12) 
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or 



and the 5-wave and D-wave wavefunctions oscillate in phase, with a ratio of 
v2. If there are nodes in the wavefunctions, then they must overlap. It is 
remarkable that such detailed information about nuclear wavefunctions can be 
deduced directly from QCD. Notice that the 1/N C expansion is equivalent to 
the short-distance expansion. Eq. (|l3|) is a remarkably accurate representa- 
tion of the short-distance part of deuteron wavefunctions in conventional NN 
models with bnth hard- and soft-cores whose long-range part is governed by 
pion exchangees. At large distances in the real world, the angular momentum 
barrier is not negligible. With Ma^ = 1/AjvjV) the tensor force dominates the 
angular momentum barrier for distances r <C 1/Ann, which are outside the 
validity of the EFT description. It might prove useful to assign a more realistic 
large- N c scaling to the angular momentum barrier and treat it as a perturba- 
tion using the WKB approximation, which becomes arbitrarily accurate in the 
large- N c limit. 

The phase (f>o is determined by matching logarithmic derivatives of the 
interior and exterior wavefunctions, 



V^Rcot(^^R)=l + f^t a nU^ + ^y (14) 
This equation determines the renormalization-group flow of the four-nucleon 



contact operatorEJIlj. The multi-branch structure of the right side of Eq. (14) 
is a consequence of the nonperturbative treatment of the tensor force. Treating 
R = 1/A as a Wilsonian cutoff, we see that the tensor force becomes pertur- 
bative when A/A^jv <C 1. In the large- N c limit, A/ Ann oc N% ^> 1 and the 
tensor force becomes nonperturbative at arbitrarily small momentum scales. 
This suggests an interesting bound-state spectrum in the large- N c limit. 



4 A Deuteron in the Large-7V C Limit? 

It is an intriguing mystery that QCD has a single scale, Aqqd, which is 
0(100 MeV), while characteristic nuclear binding energies are 0(1 MeV). In 
particular, the deuteron binding energy, —B, is B — 2.224575 MeV. The prob- 
lem might not actually be this severe. The tensor force in the chiral limit is 
governed by a single scale, F^. If we neglect geometrical factors and assume 
that the momenta of the bound nucleons is ~ F % , then we would expect a 
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binding energy of order F%/M ~ 10 MeV. Of course this argument breaks 
down as we vary N c , since the axial coupling grows with N c . Before consider- 
ing the bound state spectrum in the pionful EFT, we will consider a simplified 
scenario. 

At distance scales much larger than the pion Compton wavelength, the 
pion can be integrated out and the EFT simplifies considerably LLj. IijlEFT^) 
observables are determined purely by four-nucleon contact operatorsu. In the 
3 <Si — 3 Di coupled-channel there is a single momentum-independent interaction, 
Co oc N c . We will assume that this operator is dominant at low energies in the 
large- iV c limit. The 3 5i scattering length is then 

where we use a square-well regulator and the arrow indicates that renormal- 
ization has been performed and Co is a renormalized quantity. Hence a oc N%. 
If the sign of Co is positive, then there is a bound state with binding energy 

Hence, B oc 1 /N^. What do we expect if we include pions in the EFT? We have 
seen that the tensor force dominates the potential at long distances. Hence we 
can write the Schrodinger equation in coordinate space as 



\V>=E\q>, (17) 



where the N c dependence has been scaled out and the barred quantities are 
iV c -independent. We can scale out the N c dependence through a rescaling of 
the coordinate, r — > N£r. The eigenvalue equation then becomes 

\*b>=B\* b >, (18) 



iV c 5 \M f 3 , 

and we again find B oj_L/iV c 5 . This argument can be formulated in a more 
physical way as follows til We can write the potential as 

Nr 

Vt - (19) 

where tb can be loosely interpreted as the size of the bound state and I have 
assumed that all dimensional parameters are given by powers of Aqcd, which 
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is defined to be unity. Now assume that there is a bound state in the 3 Si 
channel at large- iV c . The presence of a bound state requires a balance between 
the kinetic energy, N c v B /2, and the potential energy, and therefore we expect 
Vt ~ N c v B . Here vb is the velocity of the bound state. As required for 
consistency, N c drops out of this relation and we arrive at vb 2 ~ rs~ 3 - This 
is simply a statement of the virial theorem. Since the nucleon momentum is 
N c vb, it follows from the uncertainty principle that vb ~ (iV c rs) _1 - Again 
using the virial theorem, we find 

r B otN 2 v B otl/Nl B(xl/N c , (20) 

and we have a bound state that is large, slow, and shallow in the large--/V c 
limit. Hence our naive scaling expectations in the pionful EFT agree with 
those from EFT(/f). Unfortunately both of these arguments would appear to 
be too naive as they ultimately rely on the virial theorem constraint that the 
coordinate scales as N 2 , or, equivalently, that the momentum transfer scales 
as 1/N 2 . As we will see below, a more detailed study in the pionful theory 
suggests that there are an infinite number of bound states in the large- N c limit, 
with the ground state binding energy increasing as a complicated function of 
N c . This last point is not surprising, given the large- N c wavefunctions found 
in the previous section. If one acts with the kinetic energy operator on the 



wavefunction of Eq. (12), one finds a complicated (harmonic) dependence on 



5 The Bound-State Spectrum 

One way of studying the bound-state spectrum is via the scattering length. 
Although a scattering length is strictly speaking not defined for a 1/r 3 poten- 
tials, we can impose an infrared cutoff, X, on Vt and match to a wavefunction 
at r > X of the form r — a Sl . This gives 

X (—1 + 4^6Afa 7r /A' tan (2^/QMa^/X + fo)) 

a Sl = ^ z-L, (21) 

a v /X tan a v / X + (f>Q J J 

The purpose of this formula is to give analytical understanding of the bound 
state spectrum in the large- N c limit. We would expect the infrared cutoff 
to take a value around l/m w . Choosing an ultraviolet cutoff R = 0.375 fm 
and a square-well depth of MVq = —2.07 x 10 6 MeV 2 which reproduce the 
scattering length for the physical value of N c = 3 in the full problem, gives 
0o = —7.6. I then choose the infrared cutoff X = 0.97 fm to reproduce 
the scattering length in the approximate formula. In order to explore the 
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Figure 1: The 3 Si scattering length as a function of N c . The solid line is a numerical solution 
of the Schrodinger equation with a small but finite pion mass and with R = 0.375 fm and a 
square-well depth of MVq = —2.07 X 10 6 MeV 2 tuned to the physical scattering leng th at 
N c = 3 (4>o = —7.6). The dashed line is the approximate analytic solution, Eq. (Gil) with 
<fio = —7.6 and X = 0.97 fm which reproduce the physical scattering length at N c = 3. 



bound state spectrum as a function of N c , I naively scale all quantities in the 
approximate formula and in the exact numerical result by the expected large- 
Ac scaling. The approximate formula for the scattering length and the exact 
numerical result are plotted as a function of N c in Fig. (Q). The approximate 
formula is reasonably accurate in the vicinity of N c = 3 (for several branches 
of the tangent). The corresponding plot of the binding energy is given in 
Fig. (Q). The plotted expression is derived from the approximate formula 
for the scattering length using effective range theory, and has been verified 
numerically along the plateaus using the exact solution. With the given choice 
of cutoff and square- well depth, there are no nodes in the wavefunction and the 
deuteron (at N c = 3) is, by construction, the ground state. As N c is increased 
from its physical value, the deuteron binding energy grows as a complicated 
function of N c and quickly passes beyond A x and out of the EFT. There is then 
no bound state in the EFT. As N c approaches 4, a new bound state, the first 
excited state, appears at threshold and then undergoes the same decoupling 
as N c is increased. This cycle repeats itself an infinite number of times in the 
large- N c limit. The excited levels can be studied in more detail through the 
singularities of the scattering length, using the_approximate formula, Eq. (^l|). 
One can estimate the density of eigenstates EJ for the 1/r 3 potential using 
WKB; with excitation number, n, one finds \E n \ ~ tt^Kqcd with n = 1, 2 . . .. 
Therefore as N c — > oo there are an infinite number of bound states, infinitely 
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Figure 2: The deuteron binding energy as a function of N c abstracted from the approximate 
analytical solution, Eq. (|2ll), and verified numerically along the plateaus. 



dispersed in energy. But there would appear never to be more than one bound 
state present in the EFT. 



6 Caveats 

There are many caveats to the large- _/V c picture of NN scattering presented 
in this paper. Among them is the absence of the A degrees of freedom. NN 
scattering becomes a significantly more complicated coupled channel problem 
when the As are degenerate with the nucleons 0. It would be interesting to 
repeat the analysis presented here with the complete contracted spin-flavor 
SU(4) multiplet. A more serious problem lies in the fact that the large- N c limit 
in which the momentum transfer is held fixed does not experience scattering, 
and potentially does not exist as a sensible limit of QCDEJJO. One may also 
seriously question the procedure of varying N c slightly away from its physical 
value to deduce information about the large- N c limit. It is conceivable that 
the true large- N c Schrodinger equation behaves quite differently. On a related 
note, there are several issues of noncommutativity of limits. First there is 
the issue of whether to take the large- N c limit at the level of the Schrodinger 
equation, or at the level of its solution. There is also the issue of whether the 
chiral limit commutes with the large-7V c limit; note that the tensor force in 
the chiral limit is highly singular at all distances, and this has an important 
effect on the bound state spectrum. Anthropic arguments suggest that the 
fine-tuned nature of the deuteron binding energy is essential to the emergence 
of life. The results presented in this paper suggest that the large- N c limit can 
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be relevant to nature only at discrete values of N c where there continues to 
be a bound state with deuteron quantum numbers and binding energy within 
anthropic bounds. Whether this be the ground state or the iV c th excited state 
is irrelevant in the EFT. Other (large) values of N c -away from the plateaus 
in Fig. (||)- would in all likelihood not support the intelligent life necessary to 
think about large- N c QCD. 
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